Temperature Dependence of the Spin Resistivity in Ferromagnetic Thin Films 
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The magnetic phase transition is experimentaUy known to give rise to an anomalous temperature- 
dependence of the electron resistivity in ferromagnetic crystals. Phenomenological theories based 
on the interaction between itinerant electron spins and lattice spins have been suggested to explain 
these observations. In this paper, we show by extensive Monte Carlo (MC) simulation the behavior 
of the resistivity of the spin current calculated as a function of temperature (T) from low-T ordered 
phase to high-T paramagnetic phase in a ferromagnetic film. We analyze in particular effects of 
fihn thickness, surface interactions and different kinds of impurities on the spin resistivity across 
the critical region. The origin of the resistivity peak near the phase transition is shown to stem 
from the existence of magnetic domains in the critical region. We also formulate in this paper a 
theory based on the Boltzmann's equation in the relaxation-time approximation. This equation can 
be solved using numerical data obtained by our simulations. We show that our theory is in a good 
agreement with our MC results. Comparison with experiments is discussed. 

PACS numbers: 72.25.-b ; 75.47.-m 
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I. INTRODUCTION 



The electric resistivity in magnetic metals has been 
studied for years. The main effect of spin-independent 
resistivity is unanimously attributed to phonons. As 
far as spin-dependent resistivity is concerned, we had to 
wait until de Gennes and Friedel's first explanation in 
1958i which was based on the interaction between spins 
of conduction electrons and magnetic lattice ions. Ex- 
periments have shown that the resistivity indeed depends 
on the spin orientatiouf^ii^iiii^i^ Therefore, the resistivity 
was expected to depend strongly on the spin ordering 
of the system. Experiments on various magnetic ma- 
terials have found in particular an anomalous behavior 
of the resistivity at the critical temperature where the 
system undergoes the ferromagnetic-paramagnetic phase 
transition.— i^i^i^ The problem of spin-dependent trans- 
port has been also extensively studied in magnetic thin 
films and multilayers. The so-called giant magnetoresis- 
tance (GMR) was discovered experimentally twenty years 
agOf^iS, Since then, intensive investigations, both exper- 
imentally and theoretically, have been carried ontM^ 
The so-called " spintronics" was born with spectacular 
rapid developments in relation with industrial applica- 
tions. For recent overviews, the reader is referred to 
Refs. [ll| andlH. Theoretically, in their pioneer work, de 
Gennes and Friedeli have suggested that the magnetic 
resistivity is proportional to the spin-spin correlation. In 
other words, the spin resistivity should behave as the 
magnetic susceptibility. This explained that the resis- 
tivity singularity is due to "long-range" fluctuations of 
the magnetization observed in the critical region. Craig 
et ali^ in 1967 and Fisher and Langer^'' in 1968 criti- 
cized this explanation and suggested that the shape of 
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the singularity results mainly from "short-range" inter- 
action at T > Tc where Tc is the transition temperature 
of the magnetic crystal. Fisher and Langer have shown in 
particular that the form of the resistivity cusp depends 
on the interaction range. An interesting summary was 
published in 1975 by Alexander and coworkers^^ which 
highlighted the controversial issue. To see more details 
on the magnetic resistivity, we quote an interesting recent 
publication from Kataoka.-'^^ He calculated the spin-spin 
correlation function using the mean-field approximation 
and he could analyze the effects of magnetic-field, den- 
sity of conduction electron, the interaction range, etc. 
Although many theoretical investigations have been car- 
ried out, to date very few Monte Carlo (MC) simulations 
have been performed regarding the temperature depen- 
dence of the dynamics of spins participating in the cur- 
rent. In a recent work^ we have investigated by MC 
simulations the effects of magnetic ordering on the spin 
current in magnetic multilayers. Our results are in qual- 
itative agreement with measurements.^^ 

Due to a large number of parameters which play cer- 
tainly important roles at various degrees in the behavior 
of the spin resistivity, it is difficult to treat all parameters 
at the same time. The first question is of course whether 
the explanation provided by de Gennes and Friedel can 
be used in some kinds of material and that by Fisher and 
Langer can be applied in some other kinds of material. In 
other words, we would like to know the validity of each of 
these two arguments. We will return to this point in sub- 
section IIIIDI The second question concerns the effects 
of magnetic or non-magnetic impurities on the resistiv- 
ity. Note that in 1970, Shwerer and Cuddy^ have shown 
and compared their experimental results with the differ- 
ent existing theories to understand the impurity effect 
on the magnetic resistivity. However, the interpretation 
was not clear enough at the time to understand the real 
physical mechanism lying behind. The third question 
concerns the effects of the surface on the spin resistivity 
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in thin films. These questions have motivated the present 
work. 

In this paper we use extensive MC simulation to study 
the transport of itinerant electrons traveling in a ferro- 
magnetic thin film. We use the Ising model and take 
into account various interactions between lattice spins 
and itinerant spins. We show that the magnetic resistiv- 
ity depends on the lattice magnetic ordering. We analyze 
this behavior by using a new idea: instead of calculating 
the spin-spin correlation, we calculate the distribution of 
clusters in the critical region. We show that the resis- 
tivity depends on the number and the size of clusters of 
opposite spins. We establish also a Boltzmann's equation 
which can be solved using numerical data for the cluster 
distribution obtained by our MC simulation. 

The paper is organized as follows. Section II is de- 
voted to the description of our model and the rules that 
govern its dynamics. We take into account (i) interac- 
tions between itinerant and lattice spins, (ii) interactions 
between itinerant spins themselves and (iii) interactions 
between lattice spins. We include a thermodynamic force 
due to the gradient of itinerant electron concentration, an 
applied electric field and the effect of a magnetic field. 
For impurities, we take the Rudermann-Kittel-Kasuya- 
Yoshida (RKKY) interaction between them. In section 
III, we describe our MC method and discuss the results 
we obtained. We develop in section IV a semi-numerical 
theory based on the Boltzmann's equation. Using the 
results obtained with Hoshen and Kopelman's^^ cluster- 
counting algorithm, we show an excellent agreement be- 
tween our theory and our MC data. Concluding remarks 
are given in Section V. 



II. INTERACTIONS AND DYNAMICS 

A. Interactions 

We consider in this paper a ferromagnetic thin film. 
We use the Ising model and the face-centered cubic 
(FCC) lattice with size AN^ x Ny x N^. Periodic bound- 
ary conditions (PEC) are used in the xy planes. Spins 
localized at FCC lattice sites are called "lattice spins" 
hereafter. They interact with each other through the fol- 
lowing Hamiltonian 

H/ = -J^S, -S,, (1) 

where Si is the Ising spin at lattice site i, indi- 
cates the sum over every nearest-neighbor (NN) spin pair 
{Si, Sj), J(> 0) being the NN interaction. 

In order to study the spin transport in the above sys- 
tem, we consider a flow of itinerant spins interacting with 
each other and with the lattice spins. The interaction be- 
tween itinerant spins is defined as follows. 



Hm^ -'^Ki^jSi-Sj, (2) 

where is the itinerant Ising spin at position r^, and 
denotes a sum over every spin pair (si,Sj). The 
interaction Kij depends on the distance between the two 
spins, i.e. ry = \ fi — rj\- A specific form of Ki j will be 
chosen below. The interaction between itinerant spins 
and lattice spins is given by 

where the interaction lij depends on the distance be- 
tween the itinerant spin and the lattice spin Si. For 
the sake of simplicity, we assume the same form for Kij 
and li j , namely, 

Ki^j = A'oexp(-ry) (4) 
Ii.j = /oexp(-ry) (5) 

where Kq and Iq are constants. 

B. Dynamics 

The procedure used in our simulation is described as 
follows. First we study the thermodynamic properties 
of the film alone, i.e. without itinerant spins, using Eq. 
Q . We perform MC simulations to determine quantities 
as the internal energy, the specific heat, layer magneti- 
zations, the susceptibility, ... as functions of tempera- 
ture TfSS From these physical quantities we determine 
the critical temperature Tc below which the system is in 
the ordered phase. We show in Fig. [l]the lattice magne- 
tization versus T for iV^ = 8, iV^ = A^y 20. 
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FIG. 1: Lattice magnetization versus temperature T for A^^ — 
8. Tc is ~ 9.58 in unit of J = 1. 

Once the lattice has been equilibrated at T, we inject 
A'o itinerant spins into the system. The itinerant spins 
move into the system at one end, travel in the x direction, 
escape the system at the other end to reenter again at the 
first end under the PEC. Note that the PEC are used to 
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ensure that the average density of itinerant spins remains 
constant with evolving time (stationary regime). The 
dynamics of itinerant spins is governed by the following 
interactions: 

i) an electric field E is applied in the x direction. Its 
energy is given by 

He = -eE • v, (6) 

where v is the velocity of the itinerant spin, e its charge; 

ii) a chemical potential term which depends on the con- 
centration of itinerant spins (" concentration gradient" ef- 
fect). Its form is given by 

He = Dn{v), (7) 

where n(r) is the concentration of itinerant spins in a 
sphere of radius D2 centered at r. D is a constant taken 
equal to Kq for simplicity; 

iii) interactions between a given itinerant spin and lat- 
tice spins inside a sphere of radius Di (Eq. [3]) ; 

iv) interactions between a given itinerant spin and 
other itinerant spins inside a sphere of radius D2 (Eq.[2]). 

Let us consider the case without an applied magnetic 
field. The simulation is carried out as follows: at a given 
T we calculate the energy of an itinerant spin by taking 
into account all the interactions described above. Then 
we tentatively move the spin under consideration to a 
new position with a step of length uq in an arbitrary di- 
rection. Note that this move is immediately rejected if 
the new position is inside a sphere of radius tq centered 
at a lattice spin or an itinerant spin. This excluded space 
emulates the Pauli exclusion principle in the one hand, 
and the interaction with lattice phonons on the other 
hand. If the new position does not lie in a forbidden 
region of space, then the move is accepted with a proba- 
bility given by the standard Metropolis algorithmi^ 

To study the case with impurities, we replace randomly 
a number of lattice spins S by impurity spins a. The 
impurities interact with each other via the RKKY inter- 
action as follows 

= - ^ L{r,, rj)a^ ■ a-j (8) 

where 

L{ri,rj) = LoCOs(2fcF|r., - rj\)/\ri - r^p (9) 

Lq being a constant and kp the Fermi wave number of the 
lattice. The impurity spins also interact with NN lattice 
spins. However, to reduce the number of parameters, we 
take this interaction equal to J = 1 as that between NN 
lattice spins (see Eq. [T]) with however a ^ S. We will 
consider two cases a — 2 and cr = in this paper. 

III. MONTE CARLO RESULTS 

We let A^o itinerant spins travel through the system 
several thousands times until a steady state is reached. 



The parameters we use in most calculations, except oth- 
erwise stated (for example, in subsection llll BI for N^) are 
s = S* = 1 and ^ Ny ^ 20 and = 8. Other param- 
eters are Di = D2 = 1 (in unit of the FCC cell length), 
Kq =^ Iq ^ 2, Lq = 17, TVo = 8 X 20^ (namely one itin- 
erant spin per FCC unit cell), vq = 1, kp = (— )(-^)^/'^, 

ro = 0.05. At each T the equilibration time for the lattice 
spins lies around 10^ MC steps per spin and we compute 
statistical averages over 10^ MC steps per spin. Taking 
J — 1, we find that Tc ~ 9.58 for the critical temperature 
of the lattice spins (see Fig. [T]). 
We define the resistivity p as 

P=-, (10) 
n 

where n is the number of itinerant spins crossing a unit 
area perpendicular to the x direction per unit of time. 



A. Effect of thickness and effect of magnetic field. 

We show in Figs. [2] and [3] the simulation results 
for different thicknesses. In all cases, the resistivity p 
is very small at low T, undergoes a huge peak in the 
ferromagnetic-paramagnetic transition region, decreases 
slowly at high T. 

We point out that the peak position of the resistivity 
follows the variation of critical temperature with chang- 
ing thickness (see Fig. [5]) and p at T > Tc becomes larger 
when the thickness decreases. This is due to the fact that 
surface effects tend to slow down itinerant spins. We re- 
turn to this point in the next subsection. 

The temperature of resistivity's peak at a given thick- 
ness is always slightly higher than the corresponding Tc- 

Let us discuss the temperature dependence of p shown 
in Fig. 1 

i) First, p is very low in the ordered phase. We can 
explain this by the following argument: below the tran- 
sition temperature, there exists a single large cluster of 
lattice spins with some isolated "defects" (i. e. clusters 
of antiparallel spins), so that any itinerant spin having 
the parallel orientation goes through the lattice without 
hindrance. The resistance is thus very small but it in- 
creases as the number and the size of "defect" clusters 
increase with increasing temperature. 

ii) Second, p exhibits a cusp at the transition temper- 
ature. We present here three interpretations of the exis- 
tence of this cusp. Note that these different pictures are 
not contradictory with each other. They are just three 
different manners to express the same physical mecha- 
nism. The first picture consists in saying that the cusp 
is due to the critical fluctuations in the phase transition 
region. We know from the theory of critical phenom- 
ena that there is a critical region around the transition 
temperature Tc- In this region, the mean-field theory 
should take into account critical fluctuations. The width 
of this region is given by the Ginzburg criterion. The 
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limit of this "Ginzburg" region could tally with the re- 
sistivity's peak and Ginzburg temperature pi^ The second 
picture is due to Fisher-Langer— and Kataokaii^ who sug- 
gested that the form of peak is due mainly to short-range 
spin-spin correlation. These short-range fluctuations are 
known to exist in the critical region around the critical 
point. The third picture comes from our MC simulationii 
which showed that the resistivity's peak is due to the for- 
mation of antiparallel-spin clusters of sizes of a few lattice 
cells which are known to exist when one enters the crit- 
ical region. Note in addition that the cluster size is now 
comparable with the radius Di of the interaction sphere, 
which in turn reduces the height of potential energy bar- 
riers. We have checked this interpretation by first creat- 
ing an artificial structure of alternate clusters of opposite 
spins and then injecting itinerant spins into the system. 
We observed that itinerant spins do advance indeed more 
slowly than in the completely disordered phase (high-T 
paramagnetic phase). We have next calculated directly 
the cluster-size distribution as a function of T using the 
Hoshen-Kopelman's algorithmJ^ The result confirms the 
effect of clusters on the spin conductivity. The reader is 
referred to our previous work-- for results of a multilayer 
case. We will show in the next section a cluster distribu- 
tion for the film studied in this paper. 

iii) Third, p is large in the paramagnetic phase and de- 
creases with an increasing temperature. Above Tc in the 
paramagnetic phase, the spins become more disordered 
as T increases: small clusters will be broken into single 
disordered spins, so that there is no more energy barrier 
between successive positions of itinerant spins on their 
trajectory. The resistance, though high, is decreasing 
with increasing T and saturated as T — > oo. 

iv) Let us touch upon the effects of varying Di and 
D2 at a low temperatures, p is very small at small Di 
{Di < 0.8): this can be explained by the fact that for 
such small Di, itinerant spins do not "see" lattice spins 
in their interaction sphere so they move almost in an 
empty space. The effect of D2 is on the other hand qual- 
itatively very different from that of Di: p is very small 
at small D2 but it increases to very high value at large 
D2- We conclude that both Di and D2 dominate p at 
their small values. However, at large values, only D2 has 
a strong effect on p. This effect comes naturally from the 
criterion on the itinerant spins concentration used in the 
moving procedure. Also, we have studied the effect of the 
electric field E both above and below Tc- The low-field 
spin current verifies the Ohm regime. These effects have 
been also observed in magnetic multilayer ?ii The reader 
is referred to that work for a detailed presentation of 
these points. 

Let us show now the effect of magnetic field on p. As 
it is well known, when a magnetic field is applied on a 
ferromagnet, the phase transition is suppressed because 
the magnetization will never tend to zero. Critical fiuctu- 
ations are reduced, the number of clusters of antiparallel 
spins diminishes. As a consequence, we expect that the 
peak of the resistivity will be reduced and disappears at 



high fields. This is what we observed in simulations. We 
show results of p for several fields in Fig. [31 
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FIG. 2: Lattice magnetization versus temperature T with 
different thicknesses (N^) of the film: crosses, void circles 
and black triangles indicate data for Nz—5, 8 and bulk, 
respectively. Tc{bulk) ~ 9.79, Tc(iV^ = 8) ~ 9.58 and 
Tc{N^ = 5) ~ 9.47. 



B. Effect of surface 

The picture suggested above on the physical mecha- 
nism causing the variation of the resistivity helps to un- 
derstand the surface effect shown here. Since the surface 
spins suffer more fluctuations due to the lack of neigh- 
bors, we expect that surface lattice spins will scatter more 
strongly itinerant spins than the interior lattice spins. 
The resistivity therefore should be larger near the sur- 
face. This is indeed what we observed. The effect how- 
ever is very small in the case where only a single surface 
layer is perturbed. To enhance the surface effect, we have 
perturbed a number of layers near the surface: we consid- 
ered a sandwich of three films: the middle film of 4 layers 
is placed between two surface films of 5 layers each. The 
in-plane interaction between spins of the surface films is 
taken to be Js and that of the middle film is J. When 
Js = J one has one homogeneous 14-layer film. We have 
simulated the two cases where Js = J and Js = 0.2 J for 
sorting out the surface effect. In the absence of itiner- 
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FIG. 3: Resistivity p in arbitrary unit versus temperature 
T for different film thicknesses. Crosses (a), void circles (b) 
and black triangles (c) indicate data for bulk, Nz=8 and 5, 
respectively. 
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FIG. 4: Resistivity p in arbitrary unit versus temperature T, 
for different magnetic fields. Void circles, stars, black rectan- 
gles, void triangles and black triangles indicate, respectively, 
data for (a) B = 0,(b) B = O.f J,(c) B = 0.3J,(d) B = IJ 
and (e) B = 2J. 



ant spins, the lattice spins undergo a single phase tran- 
sition at Tc ~ 9.75 for Jg = J, and two transitions when 
Js = 0.2J: the first transition occurs at Ti ~ 4.20 for 
"surface" films and the second at r2 — 9.60 for "middle" 
film. This is seen in Fig. [5] where the magnetization of 
the surface films drops at Ti and the magnetization of 
the middle film remains up to T2. The susceptibility has 
two peaks in the case Jg = 0.2 J. The resistivity of this 
case is shown in Fig. [6l at T < Ti the whole system is 
ordered, p is therefore small. When Ti < T < T2 the 
surface spins are disordered while the middle film is still 
ordered: itinerant spins encounter strong scattering in 
the two surface films, they "escape", after multiple colli- 
sions, to the middle film. This explains the peak of the 
surface resistivity at Ti . Note that already far below Ti , 
a number of surface itinerant spins begin to escape to the 
middle film, making the resistivity of the middle film to 
decrease with increasing T below Ti up to T2, as seen in 
Fig. [6l Note that there is a small shoulder of the total 
resistivity at Ti. In addition, in the range of tempera- 
tures between Ti and T2 the spins travel almost in the 
middle film with a large density resulting in a very low 
resistivity of the middle film. For T > T2, itinerant spins 
flow in every part of the system. 



C. Effect of impurity 



In this subsection, we take back iV^: = Ny = 20 and 



1. Magnetic impurities 

To treat the case with impurities, we replace randomly 
a number of lattice spins S by impurity spins a — 2. We 
suppose an RKKY interaction between impurity spins 
(see Eq. [5]) . Figure [7] shows the lattice magnetization 
for several impurity concentrations. We see that crit- 
ical temperature Tc increases with magnetic impurity's 
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FIG. 5: Upper figure: Magnetization versus T in the case 
where the system is made of three films: the first and the third 
have 5 layers with a weaker interaction Js, while the middle 
has 4 layers with interaction J = 1. We take Js = 0.2 J. Black 
triangles: magnetization of the surface films, stars: magne- 
tization of the middle film, void circles: total magnetization. 
Lower figure: Susceptibility versus T of the same system as in 
the upper figure. Black triangles: susceptibility of the surface 
films, stars: susceptibility of the middle films, void circles: 
total susceptibility. See text for comments. 
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FIG. 6: Resistivity p in arbitrary unit versus T of the system 
described in the previous figure's caption. Black triangles: 
resistivity of the surface films, void circles: resistivity of the 
middle film, black squares: total resistivity. See text for com- 
ments. 



concentration. We understand that large-spin impurities 
must reinforce the magnetic order. 

In Figs, m [9] and [101 we compare a system with- 
out impurity to systems with respectively 1 and 2 and 
5 percents of impurities. The temperature of resistiv- 
ity's peak is a little higher than the critical temperature 
and we see that the peak height increases with increas- 
ing impurity concentration (see Fig. [T0|) . This is easily 
explained by the fact that when large-spin impurities are 
introduced into the system, additional magnetic clusters 
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FIG. 7: Lattice magnetization versus temperature T with dif- 
ferent concentrations of magnetic impurities Ci. Void circles, 
crosses, stars and black diamonds indicate, respectively, data 
for (a) Ci = 0%, (b) Ci = 1%, (c) Ci = 2% and (d) Ci = 5%. 



FIG. 10: Resistivity p in arbitrary unit versus temperature 
T. Two cases are shown: without (a) and with 5% of mag- 
netic impurities (b) (void circles and crosses, respectively). 
Our result using the Boltzmann's equation is shown by the 
continuous curves (see sect. IV): thin and thick lines are for 
(a) and (b), respectively. Tc ~ 10.21 for 5%. 




o I ' ' ' ' ' ' ' ' ' 1 

O 2 4 6 8 10 12 14 16 18 20 



FIG. 8: Resistivity p in arbitrary unit versus temperature 
T. Two cases are shown: without (a) and with 1% (b) of 
magnetic impurities (void circles and crosses, respectively). 
Our result using the Boltzmann's equation is shown by the 
continuous curves (see sect. IV): thin and thick lines are for 
(a) and (b), respectively. Note that Tc ~ 9.68 for Ci = 1%. 



around these impurities are created in both ferromag- 
netic and paramagnetic phases. They enhance therefore 
P- 




o I ' ' ' ' ' ' ' ' ' 1 

O 2 4 6 8 10 12 14 16 18 20 

T 

FIG. 9: Resistivity p in arbitrary unit versus temperature 
T. Two cases are shown: without and with 2% of magnetic 
impurities (void circles and crosses, respectively). Our result 
using the Boltzmann's equation is shown by the continuous 
curves (see sect. IV): thin and thick lines are for (a) and (b), 
respectively, n ~ 9.63 for 2% . 



2. Non-magnetic impurities 

For the case with non-magnetic impurities, we replace 
randomly a number of lattice spins S by zero-spin im- 
purities (7 = 0. Figures [TTl [T^] and [T3] show, respec- 
tively, the lattice magnetizations and the resistivities for 
non-magnetic impurity concentrations 1% and 5%. We 
observe that non-magnetic impurities reduce the criti- 
cal temperature and the temperature of the resistivity's 
peak. This can be explained by the fact that the now 
"dilute" lattice spins has a lower critical temperature so 
that the scattering of itinerant spins by lattice-spin clus- 
ters should take place at lower temperatures. 
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FIG. 11: Lattice magnetization versus temperature T for dif- 
ferent concentrations of non-magnetic impurities: (a) Ci — 
0% (void circles), (b) Ci = 1% (crosses) and (c) Ci = 5% 
(black diamonds). 



D. Discussion 

De Gennes and Friedeli have shown that the resistivity 
p is related to the spin correlation < Si • Sj > . They have 
suggested therefore that p behaves as the magnetic sus- 
ceptibility X- However, unlike the susceptibility which di- 
verges at the transition, the resistivity observed in many 
experiments goes through a finite maximum, i. e. a cusp. 
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FIG. 12: Resistivity p in arbitrary unit versus temperature T 
for two cases: without (a) and with 1% of non-magnetic im- 
purities (b) (void circles and crosses, respectively). Our result 
using the Boltzmann's equation is shown by the contiimous 
curves (see sect. IV): thin and thick lines are for (a) and (b), 
respectively. 




periments). We think that all systems are not the same 
because of the difference in interactions, so one should 
not discard a priori one of these two scenarios. 

In this paper, we suggest another picture to explain the 
cusp; when Tc is approached, large clusters of up (resp. 
down) spins are formed in the critical region above Tc. 
As a result, the resistance is much larger than in the 
ordered phase: itinerant electrons have to steer around 
large clusters of opposite spins in order to go through the 
entire lattice. Thermal fluctuations are not large enough 
to allow the itinerant spin to overcome the energy barrier 
created by the opposite orientation of the clusters in this 
temperature region. Of course, far above Tc, most clus- 
ters have a small size, the resistivity is still quite large 
with respect to the low-T phase. However, p decreases 
as T is increased because thermal fluctuations are more 
and more stronger to help the itinerant spin to overcome 
energy barriers. 

What we have found here is a peak of p, not a peak of 
dp/dT. So, our resistivity behaves as the susceptibility 
although the peak observed here is not sharp and no di- 
vergence is observed. We believe however that, similar to 
commonly known disordered systems, the susceptibility 
peak is broadened more or less because of the disorder. 
The disorder in the system studied here is due the lack 
of periodicity in the positions of moving itinerant spins. 



FIG. 13: Resistivity p in arbitrary unit versus temperature T 
for two cases: without (a) and with 5% of non-magnetic im- 
purities (b) (void circles and crosses, respectively). Our result 
using the Boltzmann's equation is shown by the continuous 
curves (see sect. IV): thin and thick lines are for (a) and (b), 
respectively. 



without divergence. To explain this. Fisher and Langer^^ 
and then Kataokai^ have shown that the cusp is due to 
short-range correlation. This explanation is in agreement 
with many experimental data but not all (see Ref. [lH for 
review on early experiments). 

Let us recall that < E >oc ^- < Si • Sj > where the 
sum is taken over NN (or short-range) spin pairs while 
Tx {J2^ Si)^ >= J2i.j < Si • Sj > where the sum 
is performed over all spin pairs. This is the reason why 
short-range correlation yields internal energy and long- 
range correlation yields susceptibility. 

Roughly speaking, if < Si • Sj > is short-ranged, then 
p behaves as < > so that the temperature derivative 
of the resistivity, namely dp/dT, should behave as the 
specific heat with varying T. Recent experiments have 
found this behavior (see for example Ref. |5|). 

Now, if < Si • Sj > is long-ranged, then p behaves as 
the magnetic susceptibility as suggested by de Gennes 
and Friedelii In this case, p undergoes a divergence at Tc 
as X- One should have therefore dp/dT > at T < Tc 
and dp/dT < at T > Tc. In some experiments, this has 
been found in for example in magnetic semiconductors 
(Ga,Mn)As^ (see also Ref. 15 for review on older ex- 



IV. SEMI-NUMERICAL THEORY 

In this paragraph, we show a theory based on the 
Boltzmann's equation in the relaxation-time approxima- 
tion. To solve completely this equation, we shall need 
some numerical data from MC simulations for the clus- 
ter sizes as will be seen below. Using these data , we show 
that our MC result of resistivity is in a good agreement 
with this theory. 

Let us formulate now the Boltzmann's equation for 
our system. When we think about the magnetic resis- 
tivity, we think of the interaction between lattice spins 
and itinerant spins. We recognize immediately the im- 
portant role of the spin-spin correlation function in the 
determination of the mean free-path. If we inject through 
the system a flow of spins "polarized" in one direction, 
namely " up" , we can consider clusters of " down" spins 
in the lattice as "defect clusters", or as "magnetic im- 
purities" , which play the role of scattering centers. We 
therefore reduce the problem to the determination of the 
number and the size of defect clusters. For our purpose, 
we use the Boltzmann's equation with uniform electric 
field but without gradient of temperature and gradient 
of chemical potential. Wc write the equation for /, the 
Fermi-Dirac distribution function of itinerant electrons, 
as 



df df 



(11) 



where k is the wave vector, e and m the electronic charge 



8 



and mass, e the electron energy. We use the following 
relaxation-time approximation 

{-^)coU = -{ — ), fl=fk-fl, (12) 

where Tk is the relaxation time. Supposing elastic colli- 
sions, i. e. k = k' , and using the detailed balance we 
have 



^ dt (27r)3 



Wk'Mh-fl)dk', (13) 



where VL is the system volume, Wk' ,k the transition prob- 
ability between k and k' . We find with Eq. p2)) and 
Eq. (I13p the following well-known expression 



Tk 



(27r)3 
X sin ek'^dk' 



«^fe',fe(l - cos( 



(14) 



where 9 and (f> are the angles formed by k' with k, i. e. 
spherical coordinates with z axis parallel to k. 

We use now for Eq. (fT4|) the "Fermi golden rule" and 
we obtain 



(^) 



J I < k'\V\k > p(l -cos6l)sin( 



Tk h^2Trk 

xS{k' -k)k'^dk'de (15) 

We give for the potential V the following expression 
which reminds the form of the interactions (|4|-([5| 



(16) 



where ^ (T) is the size of the defect cluster and Vq a con- 
stant. We resolve Eq. with Eq. (I16p and we have the 
following expression 

1 ^ 32V,^nmk7r ^ /■ sing(l - cosg) 

(-) = ( ^3 ) J (if2+^-.)3 (17) 

where K ^ \k — k'\ is given by 

K ^\k~k' \^k[2il- cos 9)]^ Z"^, (18) 
Integrating Eq. P?|) we obtain 



Tk 



{2knf ""^^ 



x[l 



2f2 



(19) 



where Uc is the number of clusters of size ^. We integrate 
Eq. (fTO|) for the interval [0, fci?] and we obtain finally 



1 



X 



1 + (2fcf 0' 



+ ln(l + (2fcFe)2) - 1] (20) 



Note that although the relaxation time r is averaged over 
all states in the Fermi sphere, i. e. states at T = 0, 
its temperature dependence comes from the parameters 
^ and Tic- These will be numerically determined in the 
following. If we know r we can calculate the resistivity 
by the Drude expression 



m 1 

9 ' 



(21) 



Let us use now the Hoshen-Kopelman's algorithm^ ■ to 
determine the mean value of ^ and the number of the clus- 
ter's mean size, for different temperatures. The Hoshen- 
Kopelman's algorithm allows to regroup into clusters 
spins with equivalent value for T < Tc or equivalent en- 
ergy for T > Tc- Using this algorithm during our MC 
simulation at a given T, we obtain a "histogram" repre- 
senting the number of clusters as a function of the clus- 
ter size. For temperature T below T^, we call a cluster 
a group of parallel spins surrounded by opposite spins, 
and for T above Tc a cluster is a group of spins with the 
same energy. At a given T, we estimate the average size ^ 
using the histogram as follows: calling Ni the number of 
spins in the cluster and Pi the probability of the cluster 
deduced from the histogram, we have 



(22) 



In doing this we obtain ^ for the whole temperature 
range. We note that we can fit the cluster size ^ with 
the following formula 



(23) 



where v is a fitting parameter and A a constant. These 
parameters are different for T < Tc and T > Tc- Figure 
[ni and Figure [15] show the average size and the average 
number of cluster versus temperature. To simplify our 
approach we consider that the cluster's geometry is a 
sphere with radius ^. Note that due to the fact that our 
fitting was made separately for T < Tc and T > T^, no 
effort has been made for the matching at T = Tc exactly, 
but this does not affect the behavior discussed below. 

We distinguish hereafter temperatures below and 
above Tc in establishing our theory. We write 



T<Tc 



= Po(l + C„/ncC^[-l + - 



Ml + {2kF0')]), 



{2kFif 



e = ( 



^^Y"{Tc 
Idtt 

' Bin f \ r " 



2o? ^ 



(24) 
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TABLE I; Various numerical values obtained by MC simulations which are used to plot Eqs. (|24|l and (|25|) . 



Impurity 




To 




^sup 


a 


0% S=l 
1% a=2 
2% a=2 
5% a=2 
1% a=0 
5% a=0 


9.58 
9.68 
9.63 
10.2 
9.47 
9.10 


7.4443 +/- 0.066 
7.5006 +/- 0.054 
7.7103 +/- 0.049 
7.9658 +/- 0.094 
7.2866 +/- 0.062 
7.0105 +/- 0.054 


-0.9254 +/- 0.015 
-0.9253 +/- 0.017 
-0.9856 +/- 0.016 
-1.1069 +/- 0.016 
-0.9028 +/- 0.013 
-0.9381 +/- 0.015 


-0.2267 +/- 0.006 
-0.1449 +/- 0.006 
-0.1135 +/- 0.004 
-0.0747 +/- 0.002 
-0.2106 +/- 0.006 
-0.1607 +/- 0.006 


1.51875 +/- 0.07 
1.62908 +/- 0.06 
1.64786 +/- 0.05 
2.13618 +/- 0.10 
1.51766 +/- 0.06 
1.47261 +/- 0.05 



the critical region from below. In Eq. a is the 

half-width of the peak of ric shown in Fig. 1151 po is the 
resistivity at T = and poo is that at T = oo. 

We summarize in Table 1 the different results obtained 
for the cases studied by MC simulations shown above. 
Other parameters Ainf, Bmf, Ci„/, Asup, B^up and Csup 
are fitting parameters which are not of physical impor- 
tance and therefore not given here for the sake of clarity. 
Using the numerical values of Table 1 and the average 
cluster size and cluster number shown in Fig. [TJ] and 
Fig. [ini we plot Eqs. ([M)) and by continuous lines 
in Figs. [8]fT3lto compare with MC simulations shown in 
these figures. We emphasize that our theory provides a 
good "fit" for simulation results. 



T>T, 



Poo(l + CsupricS.'^i-l + 



1 



+ ln(l + (2fcf^0')]), 
^ A 

Idtt 



1 + {2kpO' 



„/3 



Bsup exp[-£)(T - Tc)] + uq. 



(25) 



In Eq. ((241) '^^ call Tq the temperature on which the 
cluster of small size are gathering to form bigger clus- 
ter. This temperature marks the limit when one enters 




10 12 14 16 18 20 




10 12 14 16 18 20 



FIG. 14: Mean size ^ of magnetic clusters versus temperature 
T for both above and below Tc- 



FIG. 15: Number of average cluster size ric for both above 
and below Tc- 



Based on those results, we can extract the resistivity 
pi corresponding only to the addition of impurities, pi 
is defined as 

Pl = Pm ~ Pstandard, (26) 

where Pstandard is the resistivity without impurities (see 
the first line of Table 1). We compare now the resistiv- 
ity Pi with experiments realized by Shwerer and Cuddy.^ 
It is important to note that the change of behavior can 
be explained if we use the correct value for v, Tc, Tq, 
etc. Figure 1161 shows pi with magnetic impurities corre- 
sponding to the Ni-Fe system , while Fig. [T71 shows pi 
in the case of non-magnetic impurities corresponding to 
the Ni-Cr case. 



We see that the form of pi of the alloys Ni-Fe(l%) and 
Ni-Fe(0.5%) experimentally observed^ can be compared 
to the curves of 1% and 2% of magnetic impurities shown 
in Fig. [TH For Ni-Cr(l%) and Ni-Cr(2%), experimental 
curves are in agreement with our results of non-magnetic 
impurities shown in Fig. 1171 

Finally, to close this section, let us show theoretically 
from the equations obtained above, the effects of the den- 
sity of itinerant spins on the resistivity. Figure [18] shows 
that, as the density ng is increased, the peak of p dimin- 
ishes. It is noted that this behavior is very similar with 
that obtained by Kataoka.^^ In our MC simulation shown 
above, we have chosen uq — 1/4. Such a weak density has 
allowed us to avoid the flip of lattice spins upon interac- 
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-5 I ' ' ' ' ' ' ' ' ' 1 

O 2 4 6 8 10 12 14 16 18 20 

T 

FIG. 16: Resistivity pi in arbitrary unit versus temperature 
T. Void circles and black triangles indicate data for 1% and 
2% magnetic impurities, respectively. 



3% non magnetic impurities 
non magnetic impurities 




O 2 4 6 8 10 12 14 16 18 20 



FIG. 17: Resistivity pi in arbitrary unit versus temperature 
T. Void circles and black triangles indicate data for 1% and 
5% non magnetic impurities, respectively. 

tion with itinerant spins. In the case of strong density, we 
expect that a number of lattice spins, when surrounded 
by a large number of itinerant spins, should flip to ac- 
commodate themselves with their moving neighbors. So 
the lattice ordering should be affected. As a consequence, 
critical fluctuations of lattice spins are more or less sup- 
pressed, so is the peak's height, just like in the case of 
an applied magnetic field. Kataoka^^ has found this in 
his calculation by taking into account the spin flipping: 
the resistivity's peak disappears then at the transition. 
It would be interesting to perform more MC simulations 
with varying uq. This is left for a future investigation. 



200 




2 4 6 8 10 12 14 16 18 20 

T 

FIG. 18: Resistivity p in arbitrary unit versus temperature T 
for different densities of itinerants spins. Curves from top to 
bottom are for no = 0.2, 0.33, 0.5, 0.7, 1 and 1.5. 



V. CONCLUSION 



We have shown in this paper results of MC simulations 
on the transport of itinerant spins interacting with local- 
ized lattice spins in a ferromagnetic FCC thin film. Vari- 
ous interactions have been taken into account. We found 
that the spin current is strongly dependent on the lattice 
spin ordering: at low T itinerant spins whose direction is 
parallel to the lattice spins yield a strong current, namely 
a small resistivity. At the ferromagnetic transition, the 
resistivity undergoes a huge peak. At higher tempera- 
tures, the lattice spins are disordered, the resistivity is 
still large but it decreases with increasing T. From the 
discussion given in subsection IIII Dl we conclude that 
the resistivity p of the model studied here behaves as 
the magnetic susceptibility with a peak at Tc- dp/dT, 
differential resistivity is thus negative for T > Tc. The 
peak of the resistivity obtained here is in agreement with 
experiments on magnetic semiconductors (Ga,Mn)As for 
example.- Of course, to compare the peak's shape exper- 
imentally obtained for each material, we need to refine 
our model parameters for each of them. This was not 
the purpose of the present paper. Instead, we were look- 
ing for generic effects to show physical mechanisms lying 
behind the temperature dependence of the spin resistiv- 
ity. In this spirit, we note that early theories have re- 
lated the origin of the peak to the spin-spin correlation, 
while our interpretation here is based on the existence 
of defect clusters formed in the critical region. This in- 
terpretation has been verified by calculating the number 
and the size of clusters as a function of T by the use 
of Hoshen-Kopelman's algorithm. We have formulated 
a theory based on the Boltzmann's equation. We solved 
this equation using numerical data obtained for the num- 
ber and the size of average cluster at each T. The results 
on the resistivity are in a good agreement with MC re- 
sults. 



Finally, let us conclude by saying that the clear physi- 
cal picture we provide in this paper for the understanding 
of the behavior of the resistivity in a single ferromagnetic 
film will help to understand properties of resistivity in 
more complicated systems such as antiferromagnets, non- 
Ising spin systems, frustrated spin systems, disordered 
media, ... where much has to be done. 
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